QUANTUM MARKOVIAN KINETIC EQUATION 
FOR HARMONIC OSCILLATOR 



Boris V. Bondarev 

Moscow Aviation Institute, Volokolamsk road, 4-, 125871, Moscow, Russia 
E-mail: bondarev.b@mail.ru 

Specific nonequilibrium states of the quantum harmonic oscillator described by the Lindblad 
equation have been hereby suggested. This equation makes it possible to determine time-varying 
effects produced by statistical operator or statistical matrix. Thus, respective representation-varied 
equilibrium statistical matrixes have been found. Specific mean value equations have been found and 
their equilibrium solutions have been obtained. 
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1. Lindblad equation 

Statistical operator g or statistical matrix is basically applied as the quantum mechanics tool, any 
information of the nonequilibrium process proceeding within the tested system may be gained from 
[1-5]. When the process concerned proceeds within the system which fails interacting with its environment, 
statistical operator g will satisfy Liouville-von Neumann equation as follows: 

ih'g=[H,g}. (1.1) 

With provision for the fact that the system interacts with any environment, a new equation shall be 
produced [4 - 16]. Lindblad is the first one who offered the equation describing interaction of the 
system with a thermostat [10]. This work is devoted to Markovian equation, which hereby describes 
nonequilibrium quantum harmonic oscillator performance 

We will write the kinetic equation for a quantum harmonic oscillator as follows: 

ih g — [H, q~\ + ih A ^ [ a g, a + ] + [a, g a + ]^j + ih B ^ [ a + g, a] + [ a + , ga] J , (1-2) 

where 

H = hui\a + a+^\ , (1.3) 
A and B are constants. Operator a is formulated as follows: 

1 I ip 



kx . (1.4) 



Equation (1.2) is very precise to describe time-varying state of the thermostat-interacted quantum 
harmonic oscillator and its equilibrium state. 

2. Energy representation 

Now, we will define the wave functions describing specific energy state ip n (x). The very functions 
satisfy the equation as follows: 

H ip n (x) = E n ip n (x) , (2.1) 

where 

E n = fiLo(n+^\ , n = 0,l,2,... (2.2) 
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As referred to energy representation, the matrix elements of statistical operator g will be formulated 
by the equation as follows: 

Qnn' = J ¥> n ( X ) Q ¥>n' {x) dx . (2.3) 

Wave functions satisfy the following equations 



aip n = y/n(p n -! , a + (p n = Vn + 1 ip n+1 . (2.4) 

With provision for the above formulas the following matrix- formed equation (1.2) is derived: 

g nn , = -iu(n-ri) g nn > + A (2 y/(n + l)(n' + 1) Q n+ i, n '+i - (n + n') g nn ^j + 

+ B (2 vW g n _ h -(n + n' + 2) g nn ^j . (2.5) 

Now, we will writs the equation for diagonal elements of statistical matrix Q nn = W n , where W n is 
the probability referred to oscillator state ip n . The equation produced has the form as follows: 

W n = 2A((n + 1) W n+1 - n W n ) + 2 B (n W„_i - (n + 1) W n ) . (2.6) 

This kinetic equation describes particular harmonic oscillator state transitions. In this case, there may 
be gained coefficients A and B as follows 

A= ip e 5' 3 ^, B = l;Pe~ 1 i fshul , (2.7) 

where P is probability of transition per unit time; /3 = 1/kT is reciprocal temperature. 

Equation (2.6) has specific oscillator state equilibrium distribution, which satisfies the following 
equation 

A ((n + 1) W n+ i -nW n )+B(n W n ^ - (n + 1) W n ) = . (2.8) 
This equation is solved by the method as follows 

W n = (1 - q) q n (2.9) 

under the following condition 

q = — = cxp(- (3 Hlj) . (2.10) 

3. Mean value of coordinate 

Mean value b assigned by operator b is defined as 

b = Tr(bg). (3.1) 

For gaining mean value a the respective equation may be derived from formula (1.2). Using the 
equality of: 

aa + - a + a = 1 , (3.2) 

we will get the equation as follows: 

~k= ( — ioj — A + B)a . (3.3) 
Now, we can find the derivatives from mean values x and p. By applying formula (1.4) we will get: 

' P + \[kx~ = (- ilu - A + B) ( -^£= + y/nx ] . 



'm \ \Jm 

Then, we will try to equate both the real and imaginary parts of this equation: 

x = —(A — B)x+p/m. 

p = — kx— (A — B)p. 



(3.4) 



If we eliminate p from this set of equations, we can obtain the mean coordinate equation 

W+2(A- B)1c + [u 2 + {A- B) 2 ^jx = Q. (3.5) 

The above equation (3.5) provides the following solution: 

x(t) = (Ci cosLut + C 2 smujt)e- {A - B)t , (3.6) 
where C\ and C2 are arbitrary constants. 

4. Mean oscillator energy 



Now, we will find the time derivative from a+a. By applying the above equality (3.2) we will produce 
the following derivative from equation (1.2): 



a+a + 2(A- B)a+a = 2B . (4.1) 
We can define harmonic oscillator time- varying energy effects inserting the following formula in equation 

— H 1 

a+a = . 

hu 2 

Thus, the following differential equation is derived: 

77 +2(A- B)H = huj(A + B) . (4.2) 
The solution is given by the equation as follows: 

*( t ) = Gc- a <*-*)« + ^4±!, (4.3) 

where C is an arbitrary constant. 

The equation (4.2) has specific stationary solution: 

77=^4±^. (4.4) 
2 A-B v ' 

Since constants A and B are related (2.7), the stationary solution obeys the formula as follows: 

— Huj e^ hu + 1 hio , /3hoj 

H = — at. t = — T~ cth — - — . 4.5 

If it is assumed that T tents to zero, than H = Huj/ 2. If it assumed that T increases to infinity, than 
H = kT. 

5. Kinetic equation expressed in terms of coordinate and momentum operators 

We well express the equation (1.2) in terms of operators x and p. For this purpose, we will firstly 
write the equation (1.2) as follows: 

ih g — H g — q H + ih A (2 aga + — a + ag — ga + a) + ifi B [2 a + ga — aa + g — gaa + ^j . (5.1) 

Since the energy operator is equal to: 

we will insert it in equation (5.1) along with formula (1.4) to obtain the following one: 

. ■ f f> 2 KX 2 \ ^ , ( fj 2 KX 2 



2m 2 \2m 2 



i(A + B) f I , 9 ^ „ 2 \ i ~i - n ?\A 

— [p g-2pgp + gp )+k(x q-2xqx + qx ))- 

2 cj y ui ' ' 'J 

— (A — B) (x gp — p qx + ifi q) . (5-3) 

6. Coordinate representation 

As referred to coordinate function, the statistical matrix is represented by formula g(t, x, x'). In 
this case, there shall be produced the following coordinate and momentum operators: 

x = x , p = — ihd x . 

Using the above values we can write equation (5.3) by the formula as follows: 

+ (A- B)(l + xd x , +x'd x ) g. (6.1) 
As concerns statistical matrix physical interpretation, the following formula 

w(t, x) — g(t, x, x) (6-2) 

is applied for getting specific probability coefficient. 
Now, we will add new variables 

xi = — (x + x') , x 2 = x — x' . (6.3) 

In this case ^ 

d x = — di + d 2 , d x > = — 8-l - d 2 ■ 

Referring to statistical matrix g(t, x\, x 2 ) and using the above new variables we will gain the 
equation as follows: 

ih „ „ in A + B ( h 2 



d t g= — di d 2 g r- x\ x 2 g H - — d 1 g - — x 2 g + e ( 1 + x\ d\ - x 2 d 2 ) g , (6.4) 

m n fiuj \2m 2 / 

where 

In this case 

g(t, x, 0) = w(t, x) . (6.6) 
We will try solution of equation (6.4) as follows: 

g(t, X\, x 2 ) = — — f f(t,k,x 2 )e ik ^dk. (6.7) 
2ir J 

Reciprocal transformation 

fit, k,x 2 ) = J g(t, Xl , x 2 )e- lk ^ dx x . (6.8) 
With provision for equation (6.6) we will obtain: 

f(t, 0, 0) = J g(t, x, 0)dx = J w(t, x)dx = l. (6.9) 

Thus, in view of function (6.8) the following equation is formed: 

dtf ^^hl dx f + ^ xdk f-A±±(^L + n^ + £{kdk + xdx) \ I . (6 .io) 

m n aw \ 2m 2 / 
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This equation has an equilibrium solution which satisfies the both formulas as follows: 

- — d x f+£xd k f = 0, (6.11) 
m n 

\^f+ !1 ^-f + e{kd k +xd x )f = Q. (6.12) 
2m 2 

We will write the performance equation of the above formula (6.11): 

m dx dk 



h 2 k kx 



This equation has the solution as follows: 



h 2 k 2 — 2 



K X 

H — = const . 



2m 2 

This formula implies that the general solution of equation (6.11) takes the form as follows: 

/ = f(E) , 

where 

h 2 k 2 nx 2 

E = 1 . 

2m 2 

Now, we will insert this function in equation (6.12) to gain the following formula: 

df f 

— + — = . 

dE 2e 

With provision for condition (6.9) this equation has the following solution: 

f{E)=exp(-^- 



2e, 

Thus, the equilibrium solution of equation (6.10) takes the form as follows: 

/(*, = exp^- — ^_ + — J J ■ (6-13) 
We will find the equilibrium statistical matrix by formula (6.7) to obtain the following equation: 



Q(xi, x 2 ) = J ° xp (^ 



2e V 2m 

On integrating the following formula is produced: 



1 ( H k K X n \ \ ; / . 

- ! 1 e tkXl dk . 



Q(xi, x 2 ) = J— exp - ax\ — -), (6.14) 



where 

f3huj m <jj 

a = a th , (7 = — - — 

2 h 

Using formulas (6.3) we will get the equation as follows: 



/ /x [a ( a(x + x') 2 o 2 {x-x') 2 \ 
6{x,x ) = J- expl i — ). (6.15) 
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Using formula (6.2) we will get respective equilibrium probability coefficient [17] 

w(x) = cxp^— a x 2 . (6.16) 

7. Momentum representation 

As referred to momentum function, the coordinate and momentum operators are represented by the 
formulas as follows: 

x = ihdp , P = P ■ 

In view of momentum representation, the statistical matrix is formulated as g(t, p, p'). This form may 
be applied for obtaining equation (5.3) as follows: 

-(A-B)(l+pd pl + P 'd p )g. (7.1) 
As concerns physical interpretation of statistical matrix g(t, p, p'), the following formula refers to 

w{t, p) = g(t, p, p) (7.2) 

probability density applied for detecting the state when an oscillator have impulse p. 

8. Wigner function 

In order to better appreciate the physical significance of various kinetic state summands we will 
derive the equation for Wigner function w — w(t, x, p), which is specified as a quantum analog of 
classical distribution function and may be defined by applying statistical matrix g — g(t, x, x') specified 
by the relation as follows: 

w(t,X,p) = -^ jg(t,x+^fiq,x-^fiq)e- l P q dq. (8.1) 

If the statistical matrix depends on x\ and X2, than 

w(t, x, p) = — -— [ g(t, Xl =x, x 2 )e- l -P x -/ h dx 2 . (8.2) 

Z 7T n J 

Reciprocal transformation 

g(t, Xl , x 2 ) = J w(t, xi, p)e lpX2 / h dp. (8.3) 

Since there is formula (6.9) 



g{t, xi, 0) dx\ = 1 , 
Wigner function satisfies specific normalization requirement 



/«(*,,,,)*,*-!. (8.4) 



For the purpose of Wigner function, we will derive the following formula from equation (6.4) formulated 
for statistical matrix g{t, x\, x 2 ) 



p » , a , a + b ( h 2 2 h 2 K 

— o x w + nx o v w + — — - — o„w H — 

m p 2Huj \2m x 2 



d 2 w + e(2 + xd x +pd p )w^j . (8.5) 
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The equation produced is much different from its quantum analog of Fokker-Planck equation. Summands 
containing derivatives d 2 w and d 2 w may be interpreted as those describing phase space diffusion. And 
still, it is necessary to add that it is rather hard to appreciate physical significance of the formula in 
parentheses that follows coefficient e. 

Equilibrium solution of equation (8.5) should at the same time refer to the following equations: 

— — d x w + k x d„w = , (8.6) 

TO 



i—d 2 x w+ d 2 w + e(2 + xd x +pd p )w = 0. (8.7) 

Z 777 Zi 

General solution of the equation (8.6) derives by the function as follows: 

w = w(E) , (8.8) 

where 

2m 2 ' 

We will insert the above function in equation (8.7). Thus, we will get the following equation: 

d 2 w ,„ ^ N dw 
dEP +{1 + » E) dE 

where 



E ^2+( l + ^E) — + M ^-0, (8.9) 



2 Bhio 

^=_th^-. (8.10) 
Solution of this equation derives the function as follows: 

w(E)=Ce~ liE . (8.11) 

Wigner function may be derived by applying formula (8.2), providing that specific equilibrium 
function (6.14) is inserted in. In view of this we shall obtain the following equation: 

w(x, p) = — !-r- [ exp( -ax 2 - e- ipX2 / h dx 2 . (8.12) 
2 7r n V n J \ 4 a J 

On integrating we will get the equilibrium function 

This function may be formulated by the equation as follows: 

w(x, p) = \J~ cx p(~ ax 2 ^j xjl exp^ — j p 2 ^ , (8-14) 

where 

7=-^-- (8-15) 
2 TO 

Then, we will find the mean value 

x 2 p 2 = J x 2 p 2 w(x, p) dx dp . (8.16) 



The above computation gives the following formula 

4 V eJ 3hu - 1 



n 2 ( e l 3hul + 1 \ 2 
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This formula has the following result. If T = 0, the uncertainty is equal to x 2 p 2 = h 2 /4. If T — > oo, 
than a; 2 p 2 = m (kT) 2 / k. 

9. Conclusion 

The equation proposed by Lindblad for the purpose of the statistical operator describing nonequilibrium 
state of quantum harmonic oscillator is hereby considered. Initially, the statistical matrix equation 
in energy representation and diagonal matrix element equation have been derived from the equation 
concerned. Specific formulas appreciating physical significance of Lindblad equation coefficients have 
been formulated. Then, the mean coordinate equation has been derived to find any general solution. It 
was demonstrated that the mean coordinate exponentially decays in time. The mean oscillator energy 
equation has been derived to obtain the general solution and mean equilibrium energy value has been 
found. Lindblad equation has been formulated by applying coordinate and momentum operators. The 
coordinate representation statistical matrix equation has been obtained. The equilibrium statistical matrix 
formula has been derived from this equation. The momentum representation statistical matrix equation 
has been formulated. Wigner function equation has been obtained and the respective equilibrium state 
function has been found. Various temperature uncertainty relations have been found by applying Wigner 
equilibrium function. 
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